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Abstract. I review recent studies about the gravitational potential and stellar dynam- 
ics of M87 in particular, and the dynamics of the stars in the presence of a super-massive 
central black hole, in general. 

At large radii, investigations of both the X-ray emitting gas and the velocity dis- 
tribution of globular clusters indicate the presence of large amounts of non-luminous 
matter, possibly belonging to the inner parts of the Virgo cluster. 

At small radii, there is no evidence from the stellar kinematics, at most a hint, for 
the existence of a central point mass, whereas the gas dynamics reveal the presence of 
a highly concentrated mass in the centre of M87, possibly a super-massive black hole 
(bh). Given the existence of such a central mass, the stellar kinematics indicate a strong 
tangential anisotropy of the stellar motion inside a few arcseconds. The implications 
of this result for the evolution and formation history of M87 and its central bh are 
discussed. I also discuss in more general terms the structural changes that a highly 
concentrated central mass can induce in its parent galaxy. 

1 Introduction 

According to their observable properties, elliptical galaxies can be divided into 
two classes. This dichotomy is most clearly revealed in the central brightness 
distribution (cf. Fig. 3 of Gebhardt et al. 1996). Consequently, the two classes 
are commonly called 'power-law galaxies' and 'core galaxies' (though - as so often 
in astronomy - the latter term is highly misleading: these galaxies actually do not 
have a core of constant density). Core galaxies have shallow central luminosity 
density profiles with j oc r -7 , 7 5a 1.3, ellipticities EO to E3-4, elliptic to boxy 
isophotes, and negligible rotation u rot <C a. They are on average bright {My & 
— 19.5) and often radio-loud and X-ray-active, possess extended stellar envelopes 
and rich (N > 2000), extended globular cluster (GC) systems which are multi- 
modal in their properties. These galaxies are thought to be of round to triaxial 
shape supported by anisotropy in the stellar motions. 

Power-law galaxies, on the other hand, have steep central density cusps with 
7 il —1.5, ellipticities up to E7, elliptic to disky isophotes, and significant rota- 
tion velocities v TOt ~ a. They are on average fainter {My ^ —21.5) and show no 
radio or X-ray activity; their surface density follows a de Vaucouleurs profile and 
their GC systems are poor (N 5s 1500) and with a profile following that of the 
stellar light. These galaxies are believed to be of near-oblate shape supported by 
rotational motions. 

Clearly, M87 having a shallow density cusp (7 w 1.2), round isophotes, 
neglible rotation, radio and X-ray activity, an extended stellar envelope, and 
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a rich and bi-modal globular cluster system is a generic representative of the 
class of core galaxies. It is generally believed, that the core galaxies are formed 
by one or more major merger events. In face of this hypothesis, it is important 
to ask whether M87 is consistent with being a merger remnant. 
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Fig. 1. Velocity dispersion for M87 as measured for stars (filled circles, van der Marel 
1994) and planetary nebulae (open squares, Cohen & Ryzhov 1997). The full line is the 
rotational velocity measured by the latter authors. 



2 The Matter Distribution at Large Radii 

The presence of dark matter around spiral galaxies was established by studies 
of the motions of their gaseous disk, whose emission lines can be traced to large 
radii. Elliptical galaxies in general have very little cold gas and one has tradi- 
tionally used stars as dynamical tracer population. However, at large radii the 
stellar absorption-line spectra are very hard, if not impossible, to observe with 
useful signal-to-noise ratio. Thus, in order to probe the potential at large radii 
of elliptical galaxies, one needs different kinematical tracers. Possible candidates 
are gas rings, planetary nebula; (PN), and globular clusters (GC). Gas rings are 
rather rare among elliptical galaxies, and in general one is left with PNs or GCs. 
The problem with using these or stars as tracers is that they do not move on 
circular orbits like the gas, but form dynamically hot systems, which complicates 
the interpretation of the measured kinematics. 
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Fig. 2. Mass enclosed in radius r as derived from studies of stellar kinematics (open 
circles, Sargent et al. 1978), globular clusters (triangle, Mould et al. 1990, and filled 
circles, Cohen & Ryzhov 1997), and X-ray gas (thin lines: lower and upper limits, 
Nulsen & Bohringer 1995, corrected for the difference in the adopted distance to M87). 
The solid line represents a power-law fit to the filled circles: M oc r l r . 



Recently, Cohen & Ryzhov (1997) have studied the GC system of M87. The 
derived rotation velocity w ro t and velocity dispersion a are displayed in Fig. 1 
together with stellar kinematical data for the central parts. There is a clear 
change in the kinematic properties at about 100" from the centre: in the inner 
parts v Iot ~ and a decreases from « 400km s -1 to « 200km s _1 ; in the outer 
parts a increases reaching w 400kms _1 at 400" and v Iot becomes significant. 
This kinematical behaviour is very similar to that observed for NGC 1399 (cD 
galaxy in the Fornax cluster) by Arnaboldi & Freeman (1996) using PNs as 
tracers (their Fig. 2). An obvious explanation is that the change in kinematical 
properties constitutes the transition from the highly concentrated galaxy to the 
less concentrated and dark-matter-dominated Virgo cluster. 

Cohen & Ryzhov have analysed their data in terms of the mass distribution 
and found that in the range 3kpc < r < 30kpc (using a distance to M87 of 
15Mpc) their data imply a mass density p oc r~ 13 and a mass of w 3 x 10 12 M Q 
inside 44kpc. These findings are in agreement with estimates derived from X- 
ray observations (Fabricant & Gorenstein 1983, Nulsen & Bohringer 1995), see 
Fig. 2. It is intriguing that the derived density slope of 7 = 1.3 agrees well with 
1.4 recently predicted for the inner parts of dark matter halos by high-resolution 
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simulations of CDM cosmogonies (Moore et al. 1997). 

3 Small Radii: Stellar Kinematics and the Black Hole 

Because of its stellar kinematics, M87 has for a long time been suspected to 
harbour at its centre a super-massive black hole 1 (bh), though this conclusion 
was always controversial (cf. Sargent et al. 1978, Binney & Mamon 1982, Merritt 
1987). The best currently available photometry (Lauer et al. 1992) and stellar 
kincmatical data (van der Marel 1994) show a weak density cusp and a slightly 
centrally peaked velocity dispersion a (see Fig. 1). As van der Marel's analysis 
showed, these data can easily be interpreted by a pure stellar model with constant 
mass-to-light ratio and anisotropy 2 of w 0.5 (van der Marel 1994), i.e. radially 
anisotropic as expected for a galaxy formed by violent relaxation. Alternatively, 
one can explain the centrally peaked a by a central bh of a few billion solar 
masses. However, such a model does not work with radial anisotropy near the 
bh, which would give too high a central a, and requires the opposite: isotropy 
or tangential anisotropy, depending on the mass of the bh. 

Bender et al. (1994) found a change in the profile of the stellar absorption 
lines: the coefficient h 4 of the Gauss-Hermite expansion of the profile changes 
sign at r « 3", such that h\ < inside and hi > outside. This implies a 
change in the underlying dynamical properties in the sense that f3 r <3" < Pr>3" , 
as predicted by models with central BH (though it would be very hard to quantify 
this). 

Actually, since 1994 we know from HST observations of a central gas-disk that 
M87 hosts a central bh of M. = (3.2±0.9) x 1O 9 M (Harm et al. 1994, Ford et al. 
1994, Marconi et al. 1997 and this volume). Starting from this fact and assuming 
spherical symmetry and a constant stellar mass-to-light ratio, Merritt & Oh 
(1997) solved simultaneously the projection equation for the stellar kinematics 
and the Jeans equation describing dynamical equilibrium. This gave them a non- 
parametric estimate for the anisotropy implied by the observed stellar kinematics 
and the bh. Their results for a bh of 1, 2.4, and 3.8 billion solar masses are 
displayed in Fig. 3. For a sufficiently massive bh, the central stellar motions 
are strongly tangentially anisotropic (fi ^ —1). The more massive the bh, the 

1 Clearly, from stellar dynamical and similar arguments, one cannot infer that the 
dark object is of the size of its Schwarzschild radius. However, lacking other plausible 
explanations for the high mass concentrations at the centres of many galaxies, it is 
generally believed that these are actually black holes, and I adopt this hypothesis 
throughout this paper. 

2 The anisotropy parameter is defined to be 

P = 1 - Ot/o 2 r , 

where a t and oy denote the tangential and radial velocity dispersions. For an isotropic 
system (3 = 0; for radial anisotropy /3 > reaching f3 = 1 for pure radial orbits; and 
for tangential anisotropy /3 < with /3 = — oo for pure circular orbits. 
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Fig. 3. Velocity dispersions and anisotropy vs. radius for M87 as inferred by Merritt 
& Oh (1997) from the stellar kinematics for a bh with assumed mass of (a) 1, (b) 2.4, 
and (c) 3.8 billion solar masses. Dotted lines give 95% confidence bands. 



stronger the inferred anisotropy and the larger is the radius at which /3 changes 
sign. Outside of this radius, the stellar motions are slightly radially anisotropic 
in agreement with van der Marel's (1994) study. For the most recent estimate of 
M # , the effect is very significant leading to highly tangential anisotropy. In the 
following section, I discuss possible explanations for this effect. 

4 Stellar Dynamics Around a Black Hole 

Can we understand the abrupt change in the stellar motions from radial to 
tangential anisotropy near the BH? Let us consider various possible scenarios for 
the formation of the bh at the centre of a galaxy. 

Adiabatic Growth of the Black Hole 

If the bh grows slowly over a long period of time, the actions of the stellar 
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orbits are conserved. However, the mapping between actions and ordinary phase- 
space coordinates changes as the potential evolves from harmonic (for an initially 
isothermal core) to Keplcrian. This change results in a stellar density cusp with 
p oc r~ 3 / 2 and a mildly tangential anisotropy of at most (3 = —0.3 (Young 1980, 
Goodman & Binney 1984, Quinlan et al. 1995) inside the sphere of influence of 
the BH, which has radius 

r 7l = ^(~l"forM87). (1) 

While this radius is of the correct size, the effect is much too weak in order to 
explain Mcrritt & Oh's result of (3 ^ — 1 . Furthermore, M87 has a central density 
slope 7 w 1.2 < 3/2. 

Growth by the Capture of Stars 

A BH may grow by tidally disrupting and capturing stars that happen to come 
too close to it. Such stars will predominantly have low angular momenta, while 
near-circular orbits are hardly affected. Hence, this process gives rise to (3 < 0. 
However, the distance from the BH that a star must reach before it gets destroyed 
is of the order of the BH's Schwarzschild radius, which is smaller than by 
several orders of magnitude. 

Growth by the Accretion of Black Holes 

If two galaxies, each hosting a massive BH at its centre, merge, the BHs will sink 
to the centre of the remnant and finally merge as well (Ebisuzaki ct al. 1991). 
Begelman et al. (1980) have outlined the stages of this process: 

(1) Due to dynamical friction with the background stars, the bhs sink into the 
centre and form a BH binary. 

(2) The binary looses energy and angular momentum (it "hardens") due to 
three-body interactions with passing stars. 

(3) When the separation between the bhs has become sufficiently small, gravita- 
tional radiation becomes very efficient in hardening the binary until it finally 
merges to a single BH. 

Depending on the time-scale of the whole process, repeated events of this kind are 
possible 3 , in particular for a central-cluster elliptical, such as M87, which over its 
lifetime has likely cannabalised many minor companions. The least understood 
mechanism here is (2), which is also the process most relevant to the possible 
effect on the stellar dynamics. One problem, for instance, is whether or not the 
eccentricity of the binary BH increases, which in turn is relevant for (3), since 
for a highly eccentric binary gravitational radiation can take over at higher 
energies, i.e. earlier, than for a less eccentric binary. Another problem is loss- 
cone depletion: three-body interactions that harden the binary eject stars out 

3 One BH-merger has to be finished before a third BH arrives, since more than two bhs 
cannot co-exist for long in a galaxy's centre, as all but two of them will quickly be 
ejected by sling shots. 
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of the centre, and the number of interaction candidates diminishes. Eventually, 
this may even halt the hardening before gravitational radiation can take over 
(the re-fueling of the loss cone due to two-body relaxation is much too slow). 

It is clear that process (2) will create tangential anisotropy, since, as for the 
tidal disruption, low-angular-momentum stars are more likely to interact with 
the BHs - they are, however, not eaten by the bh but ejected via a sling-shot. 
Evidently, this process reaches out to the radius r/j, where the circular motion 
around the bhs equals the mean stellar velocity dispersion. 




Fig. 4. Response of a Jaffe model to the hardening of a bh binary as computed by 
Quinlan & Hernquist (1997). The dashed lines show the initial (a) density, (b) velocity 
dispersion, and (c) anistotropy. The solid lines show these quantities when the binary 
reaches a certain hardness. The binary's components have equal mass of 0.04, 0.02, 
0.01, 0.005, and 0.0025 (leftmost) in units of the total mass in stars. Reprint from New 

Astronomy, Vol. 2, Quinlan & Hernquist, "The dynamical evolution of massive black hole binaries - 
II. Self-consistent A^-body integrations", p. 533-554, 1997 with kind permission of Elsvier Science - 
NL, Sara Burgcrhartstraat 25, 1055 KV Amsterdam, The Netherlands. 



Recently, Quinlan (1996) and Quinlan & Hernquist (1997) have studied pro- 
cesses (1) and (2) by detailed simulations. They developed an TV-body code in 
which the BH-bh interaction is computed exactly (Newtonianly; relativistic ef- 
fects are unimportant before stage 3), the BH-star interactions by a softened 
Keplerian force, and the star-star interactions via an expansion of the stellar po- 
tential in basis functions, i.e. essentially collisionless. Their code thus follows the 
dynamical evolution of the bhs and stars in a self-consistent way. In their simu- 
lations, Quinlan & Hernquist started with the stars being in spherical isotropic 
equilibrium following a Jaffe model (which has 7 = 2) and various choices for 
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the masses and initial positions of the bhs. Their general conclusions relevant 
here are as follows, (i) The total mass of the stars ejected from the centre by 
three-body interactions is about twice the mass of the bh binary, (ii) Wandering 
of the bh pair significantly increases the loss cone and mitigates the problem 
of its depletion, (iii) Inside about r h , the density profile has become much shal- 
lower, almost flat, the velocity dispersion increases in a Keplcrian fashion, and 
the stellar motions arc significantly tangcntially anisotropic with (3 « — 1 (see 
Fig. 4). 

Thus, at least qualitatively, this process can explain the presence of a shallow 
stellar cusp with significant tangential anisotropy around a central black hole. 
However, there are some points where M87 does not fit smoothly into this picture. 
For example, the radius inside which the density of M87 becomes shallow (« 
10") is clearly larger than the radius inside which (3 < (« 3"), while the 
simulations indicate that these radii should be similar. Also, the central cusp of 
M87 is not completely dug out as in Quinlan & Hcrnquist's simulations. These 
discrepancies indicate that reality is more complicated, possibly involving several 
such accretion events with "small" secondary BHs and/or dissipational processes 
(e.g. star formation) as indicated by the presence of the central gas-disk. 

5 Influence of a Central Black Hole on Larger Scales 

The processes discussed in the last section influence the structure and dynamics 
of the stellar system only in the immediate neighbourhood of the BH, i.e. inside 
r/,. A central bh, however, will influence all those stars that ever come near the 
centre. Many stars in triaxial galaxies are on box orbits, which pass arbitrarily 
close to the centre after a sufficiently long time. It was argued by Gerhard & 
Binney (1985) that scattering and re-distribution of box orbits by a central bh 
would cause at least the inner parts of a triaxial galaxy to become rounder or 
more axisymmetric. 

Using a modification of the TV-body code employed by Quinlan & Hcrnquist, 
Mcrritt & Quinlan (1998) have recently studied this process numerically. They 
created a stable triaxial equilibrium model by the simulated collapse of a non- 
equilibrium configuration, and followed the time-evolution of this model when 
at its centre a point mass was slowly grown. Fig. 5 shows, for three different 
bh masses, the time-evolution of the shape of the model. In all simulations, the 
galaxy model tends to become axisymmetric, even at the half mass radius. For 
the lightest BH with M./M g = 0.003 (M g denotes the mass of the galaxy), this 
process is still ongoing at the end of the simulation, whereas for M./M g = 0.01 it 
is nearly finished. The most interesting result, however, is that for M./M g = 0.03 
the process is almost finished at the time t grow when the BH has reached its final 
mass. Further simulation of the authors with varying i grow confirmed this result: 
for M./M g = 0.03 the shape becomes axisymmetric as soon as possible, i.e. at 
t = max{£ grow , idynM}, where td yn (r) denotes the dynamical or crossing time 
at radius r. 
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Fig. 5. Time-evolution of the axis ratios b/a and c/a (full lines) and the triaxiality 
index T = (a 2 — b 2 ) / (a 2 — c 2 ) (dashed) in the simulations of Merritt & Quinlan (1998). 
The mass M. of the bh grown in i grow = 15 is 0.003 (a), 0.01 (b), and 0.03 (c) of the 
total mass in stars. The numbers in the left corner of each frame are the fraction of 
particles, ranked by binding energy, that were used in computing the shape parameters. 



This behaviour can be explained by the BH inducing chaos in the orbital 
motions of the triaxial galaxy. For small M,, the orbits become weakly chaotic, 
i.e. they behave like a box orbit over some time, but after a sufficiently long time 
they fill their energy surfaces. For ever larger M # , the orbits become ever more 
strongly chaotic until at some critical M, the Liapunov time (the time in which 
two neighbouring trajectories diverge) equals the dynamical time, i.e. the orbits 
no longer resemble box orbits at all. In the simulations of Merritt & Quinlan 
this critical M. is about 2.5% of M g . In general, this number should depend on 
details of the triaxial configuration, but is likely to be of the same order, i.e. 
- l(T 2 M g . 

From the existence of this critical M # , Merritt & Quinlan draw a very in- 
teresting conclusion. In order for BHs to grow by gas accretion (the standard 
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model for AGNs), the gas has to reach the BH from large radii, i.e. it has to lose 
its angular momentum. In axisymmetric galaxies, the conservation of angular 
momentum along ballistic orbits renders gas- fueling of the centre very difficult. 
Thus, a bh may cut off its own gas supply by changing the shape of its host 
galaxy. If this picture actually applies, the bh mass should be no larger than 
the critical M. ~ 10~ 2 M g . A correlation in this sense is indeed observed: M. 
inferred from the dynamics of several early-type galaxies is always of this or- 
der (cf. Kormendy & Richstone 1995). In M87, M./M g is only ~ 0.5%, which 
might be a consequence of mergers that convert disks into spheroids and hence 
increasing M g (Merritt, private communications). 

Several BHs, however, are hosted by barred spiral galaxies (e.g. the Galaxy, 
NGC 1068). Tumbling bars are mainly made of stars on so-called x\ orbits, which 
avoid the very centre. Hence, the mechanism working on box orbits for triaxial 
bulges may not (or not as well) work for barred spirals. 

6 Summary 

The kinematics of M87 are well studied, which make this galaxy a good test 
case for the theories of galaxy formation. Outside ~ 100", the velocity dispersion 
profile rises indicating the presence of large amounts of non-luminous matter. 
The inferred density profile p oc r~ 13 is consistent with predictions from CDM 
cosmogony for the inner parts of dark-matter halos. 

The massive black hole (bh), detected in the very centre of M87 by gas mo- 
tions, together with the observed stellar kinematics implies a significant tangen- 
tial anisotropy of the stellar motions. Among the formation histories discussed 
for a bh in a galactic centre, only the model of accretion of other massive bhs, 
originating from the centres of cannabalized companions, can explain such a 
strong anisotropy. This scenario also predicts a shallow stellar density cusp as 
observed for M87. (Quantitatively, there are some discrepancies, which may well 
be due to over-simplification in the simulations of this process.) 

A massive bh at the centre of a triaxial galaxy renders, by the destruction 
of box orbits, the shape of its host axisymmetric. This mechanism becomes very 
fast once the BH mass reaches a critical value, which is of the order of 1% of its 
host's mass. Since the conservation of angular momentum along ballistic orbits 
in an axisymmetric galaxy obstructs gas-fueling of the centre, this process may 
pose an upper limit for the mass a BH can reach by gas-accretion. An upper limit 
of this order is indeed observed among BH masses inferred from the dynamics of 
early- type galaxies. 
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